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Abstract 

K> , The four-point perturbative contribution to the spherical partition function of the gravi- 

H I tational Yang-Lee model is evaluated numerically. An effective integration procedure is due 

to a convenient elliptic parameterization of the moduli space. At certain values of the "spec- 
tator" parameter the Liouville four-point function involves a number of "discrete terms" 
which have to be taken into account separately. The classical limit, where only discrete 
terms contribute, is also discussed. In addition, we conjecture an explicit expression for 
this partition function at the "second solvable point" where the spectator matter is in fact 
another 7W2/5 (Yang-Lee) minimal model. 
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1. Introduction 

This work is a direct continuation of the previous developments of refs. |2j and j3] where a 
kind of analytic-numeric extrapolation of the perturbative series has been applied to study 
the spherical partition function of a perturbed conformal matter immersed to the quantized 
2D gravity. In these works the scaling Yang-Lee model (complemented with a conformal 
"spectator matter") is taken as an example of the perturbed matter theory. Similar develop- 
ment for the massive free Majorana fermions (the scaling Ising model) coupled to 2D gravity 
is reported in ref. |lj. 

For the further convenience we summarize in this section few main concepts of 0]. This 
is mostly to facilitate subsequent references and introduce notations. It doesn't make the 
present report self consistent, a minimum acquaintance with the previous one |3^ is assumed. 

Gravitational Yang-Lee model (GYL) can be conventionally defined through the La- 
grangian density 

-^GYL = '^matter + >Cl + £gh (1.1) 

where £l is referred to as the Liouville component of the theory 

C^ = ^{da<py + fie"'^ (1.2) 

and £gh is the Lagrangian of the Cgh = — 26 conformal ghost BC system. As usual, (j) is the 
Liouville field, /i is the cosmological constant and 6 is a parameter related to the Liouville 
central charge Cl = 1 + QQ^ through the "background charge" Q = b + b~^. The matter 
component consists of some unperturbed "spectator" CFT with central charge Csp and the 
7VI2/5 minimal CFT model (critical Yang-Lee model with cyl = —22/5) perturbed by the 
only A^2/5 non-trivial primary field ip, i.e., the basic Yang-Lee field of dimension A = —1/5 

>Cn,atter = -^sp + CyL + \-^{x)^''^ (1.3) 

ZTT 

Parameters h in ()1.2|) and g in the interaction term of eq. ()1.3|) . are determined by Cgp and A 
through the balance equations 

Csp + Cyl -F Cl + Cgh = (1.4) 

g{Q-g)^^ = \ 

In the case of spherical geometry, which we only consider in the present study, any details 
of the spectator matter are not important except for the parameter Cgp. This component 
is added to get a formal access to the parameter }? of the model, in particular, to have a 
link with the classical limit. The interaction term contains a dimensional coupling constant 
A ~ /i^/'' where 

P = bg~' (1.5) 



The spherical partition function Z{fi, A) of the model is developed as a systematic 
perturbative series in the coupling constant A 

Z(^ = Y H^ (((^e'^'t'V)) (1 6) 

Z(/i,0) ^(27r)"n!\\^^ ^//lg ^^-"^^ 

where ((((/jc^^'^) ))tp are (integrated and normalized) ra-point functions in the Liouville 
gravity (i.e., at A = 0). 

The normalized correlation functions scale as^ 

(((^e2^<^)"»,^ = (7r/.)-"/^a„ (1.7) 

with some dimensionless numbers a„, i.e., the perturbative development is in fact a series in 
powers of the dimensionless scaling parameter \n~^^^. It is convenient to introduce also the 
unnormalized correlation functions 

In the Liouville gravity these correlation functions are evaluated as the integrals of the 
products of the Liouville and A^2/5 n-point functions 

{{w''^)\^ = {CC{x,)CC{x,)CC{x,))^^ X (1.9) 

(V9(Xi) . . . V9(x„))yl (Vg(Xi) . . . Vg{Xn))i^ (fx^ . . . (fXn 



where (. . .) j^, (. . .)yl and (. . .)^ are related to respectively ghost, AI2/5 minimal model and 
(unnormalized) Liouville correlations. We also denote Vg = exp{2g(j)). Geometrically the 
n — 3 dimensional integral in ()1.9|) is the integral over the moduli space of a sphere with n 
punctures. 

As in f3', in this paper we are mainly interested in the fixed area partition function Z^(A), 
which is related to Z{fi, A) as 

Z(/i,A)= / Z^(A)e-^^— (1.10) 

The lower limit (0) here is simply a particular prescription how to regularize the divergency 
of the integral at small A 3J. The fixed area partition function has the following scaling 
form 

ZAiX) = ZAiO)z{h) (1.11) 

where Za{0) scales as A^^l^ . The scaling function z{K) is a regular expansion 



{h) = Y,Zr.{-hf (1.12) 



z ^ 

n=0 



■^Here notations differ slightly from those of j3j, where the coefficinets a„ were dimensional and included 
the multiplier (27r)"n!. 



in the fixed area diinensionless scaling parameter 

h = \(-] (1.13) 



^^yp 



TT 



Apparently, coefficients Zn are related to the numbers G„ in ()1.8|) as 



a„r(-l -b 



-2^ 



" (27r)"n!r(rip-i - b-^ - 1) 
Sometimes it is convenient to express these coefficients directly 



T 



through the non-normalized fixed area correlation functions 



1.14) 



;i.i5) 



and the fixed area Liouville partition function 

^L - n^bT{b-^ - 1) ^^-^^^ 

It was argued in [T that z(/;,) is an entire function of h. 

The infinite volume specific energy S^i^db"^) = — A^c enters the leading A ^ oo 
asymptotic of Za{X) 

A \ — Qir/6ir 



Za{\) ~ ^oo ( - 1 exp {-8^,,{b)A) (1.18) 

This quantity scales as 

/ic = /o(&')A^ (1.19) 

Here we introduced dimensionless function foib"^) which is an important universal characteris- 
tic of GYL. Numerical study of this quantity at different values of the "spectator" parameter 
6^ is the main topic of ref. jH] and of the present article. Parameters Qm = &ir + &ir and 6ir 
in p.l8|) are fixed by the "IR central charge balance" |3j 

1 + 6QJ^ + Csp = 26 (1.20) 

In terms of the scaling function z{h) eq. ()1.18j) reads as the following asymptotic behavior 
at h -^ oo 

log z{h) = nfo{b)hP + {6 + 1/2) \ogh + 0(1) (1.21) 



where 

6 + 1/2 = p{Qb-'-QinK^) (1-22) 

All the above considerations are either kinematical or based on natural physical assump- 
tions. A less trivial observation is that in GYL (as well as in a number of other important 
gravitational models, see e.g., |lj) the asymptotic ()1.2H) holds in the whole complex plane 
away from the negative real axis and that all zeros of this entire function z{h) are real. This 
property has been established numerically in the classical limit 6^ = [2] and in the special, 
exactly solvable case 6^ = 0.4. It is then extended, as a conjecture, to the whole region of the 
parameter. Combined with the asymptotic behavior p.21|) this feature leads to an effective 
analytic-numeric algorithm, which allows to restore the scaling function to an impressive 
accuracy starting from a few first perturbative coefficients Zn in the expansion ()1.12j) (see jH| 
for more details). 

In the previous study we used these coefficients up to z^, where the numerical evaluation 
of the matter and Liouville correlation functions doesn't offer any technical difficulties. In 
the present article we develop a method of numerical integration over moduli, which allows 
to calculate the four-point function in ()1.9|1 to a precision sufficient to improve the results 
of|S]. 

For completeness we quote here the expressions for Z2 and zs used in jH] 

j^~'^/'>ib')Tib-'-lhi2gb-b') 

where 

73/2(1/5) 

5?^-l-Sl«l... (1.24) 

in the last expression is related to the basic structure constant C^p^^ in the critical Yang- 
Lee model in]. Special function Th{x) is the standard element of the Liouville field theory 
construction (see or jH] for the definitions and properties). 

2. Matter four point function 

In the four-point function the integral ()1.9|1 reduces to 

((<^e25<^)')^^ = f GYL{x,x)GL{x,x)d^x (2.1) 

where both in the Yang-Lee four point function 

GyUx, x) = (<^(0)<^(1)<^(oo)<^(x))yl (2.2) 

and the Liouville one 

G^{x, x) = (V;(0)K;(l)V;(oo) V,(x))l (2.3) 

5 



7^- 


-29/6(^2)^(^-2 _ 


1)7(2(76-62) 


mn 




7^- 


81(1 + 6-2- 

-39/^(62)1(6-2 - 


- 2^6-1) 
1)7(3^6 -62)Tj 


m 




486r(2 + 6-2 - 

7^/^(1/5) 


-3gb-^)T,i3gy 
- = 1.91131... 


mg) 





we have used the projective invariance of the integrand to put Xi, X2 and X3 to 0, 1 and 00 
respectively. 

The A^2/5 four point function reads exphcitly [6] as 

Gyl{x, x) = Ti{x)Ti{x) - h'T^{x)T^{x) (2.4) 

where k is from ()1.24|) and the bfocks Tj^^pix) are expressed through the hypergeometric 
functions 

T^{x) = x2/5(l - a;)i/^2i^i(2/5, 3/5, 6/5, x) (2.5) 

J-^(x) = x'/\l - xf'SF^{l/^, 2/5, 4/5, x) 

3. Liouville four point function 

The non-normahzed Liouville four point function ()2.3|) can be represented as an integral 
over the intermediate momentum P [7j (see also [HI from where the expressions below are 
read off). In our context it looks more convenient to start with the fixed area A Liouville 
four-point function 

/ A\ (49-Q)/fe 

G^^\x,x)=[-] ^l(x,x) (3.1) 



Try 
Function gi^{x,x) in general can be evaluated through the following integral representation 

^l(x, x) = r(4^^_i^V^ - 1) / ^'^^^^^^ (^) ^' ^^^ ^^-2) 

The prime near the integral sign denotes possible discrete terms (see below) and 



^^M^TA ^» ^» 



X (3.3) 



is the general four point conformal block with all four external dimensions Ag = g{Q — g) = 
6/5, the central charge cl and the intermediate dimension Q'^/A + P^. This function was 
introduced in j3]. In Appendix A we recapitulate some details end explicit constructions 
concerning this object. As in the previous paper [^ we use the notations 



,,2M2-2i,»\'«-*''>/' T*((>)T*(2j) 



^^-(-'<y)^'-'') nni(2,-Qi2) f^'^) 



and 



(p.^ 7r^T,(2^P)T,(-2zP)T^(2^ - Q/2) 

'^ ' Tl{b)Tl{2g - Q/2 - tP)Tl{2g - Q/2 + zP)Tl{Q/2 - iP)Tl{Q/2 + iP) ^ ' ' 



The last function enters the integrand in P in ()3.2j) and therefore a quick numerical evaluation 
is important. Expression ()3.5p admits the following integral representation j3] 

rg{P) = sinh(27rr^P) sinh(27r6P) (3.6) 

dt sin^ Pt (cosh2(Q - 2g)t - e''^* cos^ Pt) 



X exp —8 



sinh6tsinh6~^t 



(convergent a.t g > Q/4:). The integral here is convergent ii g > Q/4:, i.e., in our example, at 
b'^ > bl (where b^ = {11 — 4a/6)/5 = 0.2404 . . ., see below). At smaller values of 6^ a slightly 
more complicated expression applies 

sinh(27r6-ip) sinh(27r6P) 7^(3/2 + b-^/2 - 2gb-^) 
'''^ ' ~ 72(3/2 + 6-2/2 - 2^6-1 + z6-ip)72(3/2 - 6-2/2 - 2gb-^ - ib-^P) ^ ' ' 

[^ dt sin^ Pt (cosh2(6 - 2g)t - e'^^ cos2 Pt) 

X exp -8 / ^ — — — . , , , 

Jq t smhotsmho-^t 

This integral representation converges in the interval < ^2 < (23 — 4v^l9)/15 = 0.371 . . ., 
which complements the region of convergence of ()3.6|1 (and has essential overlap with it). 

The integral in ()3.2j) is understood literally, i.e., it goes along the real axis, only if 
62 > 62 = 0.2404. . .. At 62 = 6g two double poles of rg{P) at iP = ±(Q/2 - 2g) cross the 
integration contour and must be picked up explicitly as the discrete terms. Then, at 62 = 
bl = (4^139- 43) /25 ^ 0.1664. . . the same happens with the poles at iP = ±{Q/2-2g-b) 
and so on. In general the pair of double poles iP = ±{Q/2 — 2g — nb) with n = 0, 1, 2, . . . 
shows up in the form of the discrete term at 62 < 62 where 



2 _ 4V25n2 + 60n + 54 - lOn - 33 
^" " 5(4n2 + 4n-3) ^^'^^ 

Below we will use the notation 

P„ = t{2g + nb- Q/2) (3.9) 

In the presence of the discrete terms expression (13. 2j) reads 

iVd-l 



9l{x, x)=J2 ^n{x, X) + r(4^^-l _V2 _ 1) J '^''9iP)^P (^) ^P i^) 



(3.10) 



n=0 

where A^d is the actual number of discrete terms 

A'd = Floor [v^l + 6-4 - 146-2/5 - b-^/2 + 1/2 (3.11) 



while in the last "integral" term the ordinary integration over real P is implied. The discrete 
terms are in a sense "logarithmic" 

D„(x, x) = MnJ^P^ (x) ^p„ (x) (2 Re fpSx) + f/„) (3.12) 

where, as in j3] , we introduced the logarithmic derivative in P of the general Liouville block 

fp{x)=i^\ogJ'p{x) (3.13) 

The logarithms appearing in this derivative are due to a kind of degeneracy which occur at 
equal external dimensions and leads to the double poles in the integrand of ()3.2j) . Apparently 
this effect hides nothing conceptually new and there is no point to talk about the Liouville 
field theory as of "logarithmic CFT"^. In the four point function with different external 
dimensions there are no logarithms. The constants Mn and Un in ()3.12|) read explicitly 

^^ {Ag - h-^ - h + khf 72(%6+(fc + n-l)62) 

^ j-J- {Ag - 6-1 + (A; + n - 1)6)2 y (2^6 + W)-i'^{l + {k + 1)6^) 

and 

Un = 2vb{Ag -Q + nb)- 2vb{Ag ~Q + 2nh) + Avb{2g + nh) - 2vb{Ag + 2nh) (3.15) 

n 

- 2vb{b) - 4(n + l)61og 6 - 46Ce + 26 ^ (^(-A;6^) + V^(l + kb^)) 

k=l 

In the last expression we have, as in ref. |lj, introduced the notation 

Vb{x) = -—logTbix) (3.16) 

ax 

for the logarithmic derivative of the T;,-function. This special function can be evaluated 
through the integral representation 

convergent in the strip < Hex < Q. Outside this region one of the following relations, 
whichever more convenient, can be used to render the argument to the strip of convergence 

Vb{x + 6) - Vb{x) = 6(-2 log 6 + ipibx) + ^(1 - 6x)) (3.18) 

Vb{x + b'^) -Vbix) = b'\2logb + tlj{b'^x) +^(l-6"^x)) 



^Whatever this last term means. 



4. Elliptic modular parameter 

After all these preliminaries we can turn to the integral ()2.H1 . For the fixed area four point 
function it reads 

J4 = (^{^e^'^Y)^^^^^ = 6 I GYL{x,x)gi^{x,x)d'x (4.1) 

As in ref. Q, in eq. ()4.1|) we used the symmetry of the integral under the six element modular 
subgroup of projective group, generated by the transformations x —>■ 1/x and x —>■ 1 — x. 
This group divides the complex plane of x in 6 regions, the fundamental region G ={Rex < 
1/2; |1 — x| < 1} and its 5 images. The integral in ()4.H1 is reduced to G while the factor 6 
in front of the integral takes into account the equivalent images. 

It turns remarkably convenient to introduce the "elliptic" modular parameter through 
the standard map 

.K(l-x) ,,^, 

where 

is the complete elliptic integral of the first kind. The integral ()4.H) becomes 

X4 = evr^ f \xil-x)el{q)\'GY^iq,q)gj^iq,q)d'T (4.4) 

where F = {|r| > 1; |Rer| < 1/2} is now the standard fundamental region of the modular 

group and 

00 

Osiq) = E ^"' (4.5) 

n=— 00 

is the usual 6'-series in 

q = e*"^ (4.6) 

There are two important advantages in the form ()4.4|1 . As it has been argued in refs. jlOj . 
a general four point conformal block admits a convenient recursive representation, which 
looks particularly simple in terms of the elliptic parameter q and the so called "elliptic" four 
point block. Also, while the two blocks ()2.5p are known in closed form, the Liouville block 
(I3.3|l at general P is only computed as a power series in x or q. The recursive algorithm 
in the elliptic representation gives directly the power series in q, which is argued to have 
much better convergence then that in x. In particular, in the elliptic fundamental region 
maxF \q\ = exp(— vry^) = 0.00433 . . . while in its x image maxc \x\ = 1 near the cusps of 
the fundamental region, where the convergence of the x-series is questionable. 



Elliptic representations of the four point conformal blocks is recapitulated in the Ap- 
pendix. It is also shown there that the products of the Yang-Lee and Liouville blocks read 
in the elliptic parametrization as 

x(l - x)et{q)J'^{x)J'p{x) = (16g)^"+^'+'?.- H^{q)Hp{q)H,p{q) (4.7) 

Here and below the index a = I,(p specifies the intermediate representation in the Yang-Lee 
block, in particular A/ = and A;^ = —1/5. The corresponding elliptic blocks Ha{q) can be 
found as a series in q either through the explicit expressions ()2.5|1 



J^j{x) = (16g)9/^°[x(l - x)Y/^^el^\q)Hj{q) {U 

T^{x) = il6qY/''[x{l~x)Y/'%/\q)H^iq) 



or via the recursive relation (IA.8I). In both cases we have 



.. / X 21 2 51 4 1989 fi 111489 o 1612779 ^ /. x 

Hi(q) = 1 q^ g^ + q^ q^ + g^° - . . . 4.9 

^^^ 22^ 88^ 5456^ 1789568^ 3579136^ ^ ' 

rr r ^ . 5 2 91 4 1967 g 8211 „ 405647 ,. 
^^^^ 6^ 152^ 8816^ 70528^ 987392^ 

The "central charge deficit" part 

H,,{q) = {l-q'-q' + q'' + ...y'-^' (4.10) 

is carried out in the Appendix. Finally, the g-expansion of the Liouville elliptic block is easily 
generated through recursive relation ()A.8|1 . Although the calculation is straightforward, the 
result is somewhat cumbersome at higher orders, so that here we quote it explicitly only up 
to O(g^) (the notation Am,n is explained in ()A.2|) ) 



where p = Q/2 — g. In numerical calculations below we used the expansion of Hp{q) up to 
the order q^. In our symmetric case (all external dimensions are equal) all the elliptic blocks 
are series in g^. 

5. Modular integral 

The Liouville correlation function can be separated into the "integral part" and, sometimes, 
a number of discrete terms ()3.10p . Consequently 

Afd-l 

I4 = / ^ Alisc, n + 2^int (5-1) 

n=0 

10 



For numerical integration it is convenient to further separate each term in two parts, corre- 
sponding to two matter blocks in ()2.4|) 

T — T^^^ - ^2 7'('^) 

-^disc, n '^disc, n '^ '^disc, n 

Consider first the integral part 

£t (5.3) 



The product of the elliptic blocks 

HP{q) = H^{q)H,^{q)Hp{q) (5.4) 

is developed in a power series in q 

oo 

HM)H,^{q)Hp{q) = Y,H^^\P)q^ (5.5) 



r=0 



SO that the integrand in ()5.3|1 as a double power series in q and q. In each term the integration 
in r2 = Im r can be carried out explicitly. The result is in terms of the following function 



(16) 



2 A 1-1/2 



$(Ar.O = / rfV|16g|'^gV = , \^ / cos(7r(r - /)a;)e-"^^^('^+'-+')da; (5.6) 

Jf 7r(2A + r + /) y_i/2 

Notice, that if the integral is divergent at r2 ^ cxd, this reduction automatically takes care 
of the divergency (in the sense of analytic continuation). We arrive at the series 

(„) _ 37rRg ^ (o) 

^°* r(4(76-i-g6-i);^^'^*'^ ^^^ 

where in the last sum 

poo -^ 

4"^= / r,(P)rfP5^/j(")(P)41(P)$(p2 + gV4 + A.-l,fc,L-A:) (5.8) 

"^0 fc=0 

In our symmetric case only even L contribute. Each term in ()5.8p is suppressed by a factor 
maxF \q\ and in practice the series in L converges very fast. Below we found it sufficient 
to sum up to L = 6 to reach the 8-9 digit precision (see table H]). 

The discrete terms, if any, are treated similarly. It is again convenient to single out the 
"logarithmic" part in each of the integrals J^J^ „ 

11 



where (P„ is from eq. ()3.9|) ) 

4") = 47r(g/2 -2g- nh) I {lQqf'+'^"+^'/^-^H^pJ{q) 

Jf 



;^(°) 



(16g)^n^+A.+QV4-l^W(^) 



-i-Tp^ogHpHp 



Imrt/V (5.10) 



(only the Liouville elhptic blocks Hp{q) depend on P and therefore appear in the logarithmic 
derivative) and 

Sn = Un-AiQ/2-2g-nb)\ogl6 (5.11) 

As in the case of the integral part, the expansion in q and q induces the "level" series 






(5.12) 






L=0 

oo 

The integrals are then evaluated in the same way as ()5.7|1 

L 

Kl = 4^(<5/2 -29- nh) Y, H^:\P^)Htl,{PnW {P^ + QV^ + A^-l,k,L-k) (5.13) 

fc=0 

^2 = E (^^^41(^0 - 2H'l^liPn)) H]:\Pn)^ [Pi + QV4 + A„ - 1, fc, L - A;) 



fc=0 



where 



r' (a) / 



r(°)/ 



^;^"^(^)=^^^r(^) 



(5.14) 



A new function 



^'{A,r.l)= /"rfV|16g|^^gVlmr 



(5.15) 



(16) 



2A 



{2A + r + l)H^ 



1/2 

cos(7r(r - /)x)e-"^'^^^(2^+"+') ( 1 + ti{2A + r + /) v/T^j 

1/2 ^ ^ 



dx 



has been introduced to treat the logarithmic part. The "level" series also converges very fast 
(see next section). 

6. Numerical results 

In tabled we present the numerical values of the fourth perturbative coefficient 

X4 



za 



2A{2tiYz\^^ 



(6.1) 



12 



at different values of 6^ together with the second and third ones ()1.23|) already evaluated in 
ref. |3]. The forth column contains the preliminary estimates ^4 of Za, on the basis of 2:2, 
^3 and analytic properties of zih), as explained in 3 . Also two "exact" values are produced 
at the "solvable" points \? = 0.4 and b^ = 0.3 (see sect. 8 for details). 



52 


Z2 X 10^ 


Z3 X 10^ 


24 X 10^ 


zt'^ X 10^ 


^(cxact) ^ ^q4 


0.00 


-8.92857 


22.9899 


-31.891804 


-31.8938 




0.01 


-8.83599 


22.5977 


-31.114313 


-31.1164 




0.05 


-8.43801 


20.9364 


-27.881220 


-27.8839 




0.10 


-7.86500 


18.6240 


-23.557311 


-23.5604 




0.15 


-7.18331 


16.0092 


-18.942809 


-18.9461 




0.20 


-6.35922 


13.0616 


-14.131054 


-14.1345 




0.25 


-5.34942 


9.78831 


-9.3340267 


-9.33745 




0.30 


-4.09998 


6.28732 


-4.94949021 


-4.95242 


-4.94949020548 


0.35 


-2.55378 


2.87061 


-1.61778168 


-1.61953 




0.40 


-0.71440 


0.35675 


-0.085061507 


-0.085303 


-0.08506150735 



Table 1: Numerical values for the second, third and forth order perturbative coefficients in 
the fixed area scaling function z{h). In the forth coloumn we place the estimate of the four- 
point coefficient from the sum rules, as explained in ref. ^. Exact values, where available, 
are presented for comparison. 



At b"^ < 6q the fourth coefficient contains the contributions from the integral part and 
discrete terms in (I3.1()|l 

J2^t\n) + zt^ (6.2) 



Z4 



n=0 



In table |2 the structure of the four point integral as a sum of discrete and integral con- 
tributions is illustrated numerically. At small 6^ the integral part becomes negligible while 
more and more discrete terms appear. In order, only few first of these discrete terms really 
contribute at sufficiently small b"^. E.g., at 6^ = 0.01 in fact there are as many as 49 discrete 
terms and at 6^ = (see next section) their number is infinite. In the table we quote only 
those which count at the precision level chosen (approximately 9 decimal digits). 



Finally, our results for z^ are used in the analytic-numeric procedure described in the 
first article |3I. This allows to correct the previous numerical results for the scaling function 
z{h). In particular we improve the numerical approximations for the specific vacuum energy 
parameter /o(&^) (see eq. ()1.19j) or [3j). The new numbers /g , which take into account the 
perturbative coefficients up to z^ are presented in table El and compared with the previous 
approximations as well as with the exact values /g '^ where the last are available. 
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h' 


0.00 


0.01 


0.05 


0.10 


0.15 


0.20 


zf '^ X 10^ 


0.00 


8.6 X 10-26 


0.006087 


1.7311 


2.318 


0.799 


zf\Q) X 10^ 


-11.4450 


-10.5160 


-6.88429 


-4.187 


-6.397 


-14.93 


zf\l)xW 


-15.7357 


-15.5090 


-13.9653 


-11.51 


-14.864 




zf\2) X 10^ 


-3.94363 


-4.16729 


-4.97850 


-5.770 






zf\?>)xlQ' 


-0.662455 


-0.77641 


-1.43314 


-3.829 






zf\A)xW 


-0.092114 


-0.12415 


-0.41407 








zf\b)xW 


-0.011464 


-0.01846 


-0.13010 








zf\Q)xW 


-0.001324 


-0.00265 


-0.04666 








zf\7)xW 


-1.44668 X 10-^ 


-3.76 X 10-^ 


-0.02083 








zf\S) X 10^ 


-1.51392 X 10~^ 


-5.34 X 10-^ 


-0.01443 








zf\<d) X 10^ 


-1.52956 X lO-'^ 


-7.65 X 10-*^ 










zf\lQ) X 10^ 


-1.50085 X 10-' 


-1.11 X 10-^ 










zf \ll) X 10^ 


-1.43677 X 10-* 


-1.66 X 10-^ 










zf\l2) X 10^ 


-1.34666 X 10-9 


-2.53 X lO-s 











Table 2: Relative contributions of the integral and discrete terms in the sum fj6.2|) at different 
values of 6^. 



7. Classical limit 

In ref. j2] the classical limit of the Liouville four point function has been considered. In 
particular, the "symmetric" function with four equal dimensions a admits the following 
integral representation 



gc\{x,x) 



TT 



(2s - l)ds r2(2a + s - l)T^{2a - s)T^{s) 



r(cl) 



:F^^\a,s,y):F^^\a,s,y) (7.1) 



27rz r4(2CT)r2(2s) 

The "classical block" JF('^')((T, s, y) is expressed explicitly through the hypergeometric function 

^(^i)(a,s,i/) = /-%Fi(s,s,2s,|/) (7.2) 

while the integration contour | goes up along the imaginary axis to the left from the poles of 
T'^{2a — s) and to the right from all other singularities of the integrand. We are going to argue 
that this expression is consistent with the classical limit of the fixed area four-point function 
fl3.2|) provided we identify A^ = cr and ()7.1|) with the normalized correlation function 



lim 



g-LK.x.x) 



b^-'O Z 



(A) 
L 



gci{x,x) 



(7.3) 



14 



62 


P 


5 


f(2) 

/o 


f(3) 

/o 


/o 


/.(exact) 

h 


0.00 


0.833333 


0.111111 


0.220407 


0.218156 


0.218036 


0.2179745 


0.01 


0.831646 


0.109935 


0.220318 


0.218091 


0.217959 




0.05 


0.824462 


0.106179 


0.219719 


0.217523 


0.217353 




0.10 


0.814333 


0.103698 


0.218260 


0.215950 


0.215742 




0.15 


0.802587 


0.103880 


0.215654 


0.213053 


0.212805 




0.20 


0.788675 


0.107063 


0.211308 


0.208218 


0.207920 




0.25 


0.771700 


0.113797 


0.204224 


0.200408 


0.200042 




0.30 


0.75 


0.125 


0.192522 


0.187677 


0.187218 


0.1870437 


0.35 


0.719788 


0.142247 


0.171896 


0.165563 


0.164980 




0.40 


0.666667 


0.166667 


0.125625 


0.116905 


0.116151 


0.1158596 



Table 3: Specific energy /o determined with the use of first two, three and four perturbative 
coefficients. When available, the exact values are quoted for comparison. 

To this order, let us evaluate the integral through the residues at the infinite sequence of 
the "right" poles at s = Sn, where 

Sn = 2a + n (7.4) 

This is a rightful procedure, as the large s asymptotic of the integrand shows. Thus 



Y.D\l 



cl). 



n=0 



where 



D'^:'\x,x) = iVf ) |^(^i)(a,s„,x)|' f2Re ^log^(^^)(a,s, 
' ^ \ as 



X] 






We introduced the notations 



TT 



n-1 

n 



{2a + ky 



(4a + 2n - l){n\f -Li (4^ + ^ + A; - 1)2 



and 



U'f^ = Ai){Aa + 2n) - 2^(4(t + n - I) - A'4j{2a + n) + 2^(1 + n) 



4o- + 2n - 1 
Comparing ()7.7|) and ()7.8|) with ()3.14|) and ()3.15|) respectively, it is easy to see that 

lim6[Z„ = t/(^') 
6^0 



(7.5) 



(7.6) 



(7.7) 



(7i 



(7.9) 
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Also, it is well known (see e.g. jT^) that 

lim^pf "^ "" X ) =^(=')(a,s,x) (7.10) 

62^0 \ (J (J J 

provided a and s = P'^ + Q"^^^ are kept finite in the limit. Thus, the classical limit of the 
normalized discrete term ()3.12|) coincides with ()7.6|) 

li-^^#^ = ^i^^H^>^) (7-11) 



b2^0 z 



{A) 
L 



In the classical limit the integral in eq. ()3.10|) is saturated by the infinite sequence of discrete 
terms, the integral one vanishing. This proves ()7.3p . 

The classical block can be rendered to the form (cp. eq. ()A.16|l in the Appendix) 

In the general expression ()A.16|) only the product H^^^\q) = HspHj^ allows the classical limit. 
The classical elliptic block reads 

ff(-l)(g) = r7'=-/2-25/2(^2)^16.(g2)^^^^) ^7_^3) 

In this expression 

hs{q)= (y^) 2Fi{s,s,2s,x) (7.14) 

4s(2s-3) 2 2 s (8 s2- 14 s + 9) 4 

= 1 ^ -q H ^ -q + ... 

2s + 1 ^ 2s + 3 ^ 

while 

is the standard Dedekind product (jA.lSjl and 

oo 

eo{q^) = l-2q^ + 2q^ - 2q^^ + ...= J^ {-Tq^""' (7-16) 

n.=— oo 

the usual theta series. It is also straightforward to verify directly that H^'^^\q) is the limit 
of H^pHp as 62 ^ and P^ ^ s - Q^/A. 

For our particular application in the GYL model the classical elliptic block is evaluated 
through the explicit formula 

H('Hq)=V-'''^''iqyf\q')hs{q) (7.17) 
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while the matter elhptic blocks H^ remains the same as in eq. ()4.9|) . As usual, we are going 
to use the g-expansions. Denote 



oo 



i7„(g)i/(^')(g) = ^/^l")g^ (7.18) 



fe=0 

oo 



i7.(g)^if(^')(g) = j:e)/ 



In the classical case the integral ()4.H1 is the sum of the infinite number of discrete terms 

oo 
n=0 

where 

/(^i) = 6 / D(^')(x, x)GYL(a;, x)£x = £) - n'ji^^ (7.20) 

Jg 

Now, as in the general case, we single out the "logarithmic" integral 

jt^ = evr^iV^^i) (4") + Ml") + Ki^^) (7.21) 

where, 

4") =271 f |(16g)^"+2"+"-i/i„(g)|^ImrrfV 
Jf 

MJC^ = f |(16g)^'^+2"+"-^/i^(g)f (f/^^^) - 2 log 16) d^r (7.22) 

Ki^^ = - [ |(16g)^'^+2'^+"-i|'-^(/i„(g)/i„(g))ciV 
Jf ds 

Then, each term is obtained as a "level by level" sum in L through the double q and q- 
expansions. 

oo L 

4") =2nY,Yl ^i"^ (sn)4"2fc(s„)$'(A„ + Sn-l,k,L-k) 

L=0 fc=0 

oo L 



M(") = (f/f ) - 2 log 16) 5^ Yl 4"H^n)4"i,(5n)$(A, + 5„ - 1, A;, L - A:) (7.23) 

L=0 fc=0 
oo L 

ir(") = -2 5^ 5^ h^ (s„)/iS(.„)<l>(A, + ,„ _ 1, A;, L - A;) 

L=0 fc=0 

Every component is straightforwardly evaluated with the use of the integrals ()5.6j) and ()5.15|) . 
In practical calculations (given the required precision of 9 digits) we found it sufficient to 
sum up to L = 10 in eq. ()7.23|) . 
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Through all these calculations we arrive at the contribution of n-th discrete term to the 
fixed area perturbative coefficient 

Several first contributions are presented in table Elto manifest the convergence. The discrete 
terms sum up to the number 

oo 

^(ci) ^ J2zf\n) = -31.8918039 x 10"^ (7.25) 

n=0 

already obtained earlier [2j by means of a different numerical approach. Notice, that we had 
to take into account as many as 10 discrete terms to achieve this 9 digit precision quoted. 

8. Integrable points b^ = 0.3 and 6^ = 0.4 

In this section we comment about the two exactly solvable points b^ = 0.4 and 6^ = 0.3 in 
the family of GYL models. The solution at the "pure Yang-Lee" point 6^ = 2/5 through 
the matrix model approach has been already discussed in the previous article jSj. Here we 
recapitulate the essence very brieffy. This case of the "minimal gravity" is related to the 
ffow from the tricritical to the critical points in the generic one matrix model ^Hl- The 
corresponding scaling function Zyhixji), which is interpreted (up to an overall scale) as the 
spherical partition function of the continuous gravity, is determined explicitly as 

ZY^(x,t) = u(x,t) (8.1) 



dx"^ 
through a solution u{x,t) of the following simple algebraic equation 



x = u^ -tu (8.21 



Parameters t and x are, again up to some normalization constants, the cosmological constant 
and (/^-perturbation coupling respectively. Comparing the series expansion generated by ()8.H) 
and dHIl 

r, , X ^ ,r. ■. ( . 105 ( X \ 105 / X \2 35 / X \3 \ 

Zvl(x, t) = Zv.(0, i) (l - - (^) + - (^) - - (^) + . . .) (8.3) 

with the perturbative coefficients evaluated in the field theoretic approach, it is easy to relate 

X A/yl 



t3/2 (7r/i)3/2 

where 

7^/2(4/5) 



■4) 



., = ^-^^ = 0.0845223... (8.5) 



Then, the fixed area scahng function ()1.12|) at b^ = 0.4 reads exphcitly 

{lYLhr 



In particular 



ih) = -2V^Y1 



n=0 



n\T{n/2 - 1/2) 



Z2 = -/|l = -0.00714401 



zs = ^/yl = 0.000356752 



z^ = — l^, = -8.50615 ... X 10" 
6 



(8.6) 



•7) 



The last number and 



/o(0.4) 



3^/YLy/^ 



IT 



0.11585962159187 



(8.8) 



are quoted in tables [T] and El as the corresponding exact values for this point. 

In ref. jH] it has been argued that at b^ = 0.3 our GYL model is also a version of minimal 
gravity. It arises as a particular perturbation of the minimal 7VI3/10 model coupled to the 
Liouville gravity. Therefore there are serious reasons to believe that at this point the model 
is again exactly solvable. Unfortunately yet no exact solution in the framework of the matrix 
model approach is known. However, it is very natural to expect the existence of a closed 
analytic expression, e.g., for the scaling function ()1.12j) . We conjecture the following explicit 
form, apparently motivated by more general structures discovered in ^3] (see also jH]) 



zsMh)=T{-l/3)J2 



ih,hy 



n=0 



n\T{n/3 - 1/3) 



^.9) 



where the scale factor /eg is easily figured out through the comparison with the perturbative 
coefficients of sect.l 

7(1/3) __. ^3,^0) 



/. 



0.23254... 



"^ 21/23^5/6(3/10) 
This explicit expression gives rise to the following values of the first perturbative coefficients 



Z2 = -0.0409998231725 
Z3 = 0.00628731873887 
z^ = -0.000494949020548 
Z5 = 0.0000257778956523 



(8.11) 
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We consider the comparison of these numbers with those in the corresponding row of the 
table Has a convincing support in favor of our conjecture. Analyzing the asymptotic of ()8.9p 
one finds, in addition, that at 6^ = 0.3 



4 // 
_ ( _£I 

TT V 3 



/o = -(:^l =0.187044... 
This number is produced in table El as the corresponding "exact" value. 



3/4 



(8.12) 




Figure 1: Integration contour in the representation ()8.17|) . Dashed are the wedges where the 
integrand decreases at large \u\. 



It seems also suggestive that the conjectured fixed area scaling function ()8.9|1 follows from 
the "matrix like" algebraic equation 



for the third derivative 



u — xu = —t 
^ = ^^3/io(a:,t) 



of the "grand" partition function 

^ , ,_ xi3/3^r(4n/3-13/3)(tx-4/3)n 

^3/w[x, t) — 2_^ - 



n=0 



n\T{n/3 - 1/3) 



(8.13) 
(8.14) 

(8.15) 



xi3/3 / 81 ta;-4/3 9 ^tx-^/^f _ {tx-^l^f {tx-^'^Y {tx-^/^f , 



3 V 3640 24 



20 



24 



90 



20 



Here the cosinological constant x related to /i through 



XI 



eg 



t 



8.16) 



(7r/i)4/3 2;4/3 
Finally, let us mention a convenient integral representation for the scaling function ()8.9|) 



^3/10 (^) = 3 / exp (m'^ + /i/egM ^) u^du 
Jc 



(8.17) 



where the integration contour C goes as it is shown in figC] 
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A. Elliptic four point block 

Here we summarize the results of refs. |10j. Consider a conformal theory with central charge 
c. It will prove convenient to introduce the "Liouville like" parameterization in terms of b 



c-1 
6 



{b-' + bf 



(A.l) 



This is not a restriction for the value of c since we allow b to be complex if needed. It is also 
convenient to introduce the notation 



A. 



mb ^ + nb 



(A.2) 



Let Aj, i = 1,2,3,4 and A be the external and intermediate dimensions in the four point 
block jH], as it is illustrated in the picture below 



-^A 



Ai A3 
A9 A. 



X 



Xj(=x) 



^2(=0)/a 




Xs(= 1) 



^4 (= °^) 




(A.3) 
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According to jTU] it can be written as 

Ai A3 



-^A 



Ao A, 



X 



(A.4) 



where q is related to x as in eqs. ()4.2|) . ()4.6|) and 6*3 (g) is defined in (14 .51) . In the function 

Ai A3 



H^{q) = Ha 



Ao A4 



(A.5) 



we suppressed the dependence on the external dimensions, which are conveniently parame- 
terized in terms of Aj, z = 1, 2, 3, 4 as 



A,; 



c-1 
24 



+ A 



The i7-function (the elliptic four point block) is a power series in q 



Ha 



Ai A3 

Ao Xa 



L=0 



Ai A3 

A2 A4 



A g^ 



(A.6) 



(A.7) 



which is believed to converge at |g| < 1. It can be effectively calculated through the recursive 
relation |Tn| 



H, 



Ai A3 
A2 A4 



q 



' + ,^,A^^-"(t A4 '^- 

(m.n) 



Ai A3 

m,n+mn \ \ \ 



Here the sum is over all pairs (m, n) of positive integers and 

A — X2 _ \2 



q] {A.t 



(A.9) 



are the dimensions of degenerate representations of the Virasoro algebra with the central 
charge c. The multipliers Rm,n in (!A.8|1 read explicitly 

J, ,.. nr,s(Ai + A2 - Ar,s)(Ai - A2 - Ar,s)(A3 + A4 - Ar,s)(A3 - A4 - Ar.,s) , , 

^m,n{Xi) = I —7 (A.iUj 

llfc,«Afc,/ 

The products in (jA.lfljl are over the following sets of integers (r, s) and (/c, /) 



and 



r = —m + 1, —m + 3, . . . , m — 3, m — 1 
s = —n + 1, —n + 3, . . . , n — 3, n — 1 



k = —m + 1, —m + 2, . . . ,m — l,m 
I = —n + 1, —n + 2, . . . , n — 1, n 



(A.ll) 



(A.12) 
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while the prime sign near the last product symbol Ylk i means that the two pairs (A;, /) = (0, 0) 
and (r/i, n) are missing. 

Relation (jA.8j) leads, in particular, to a recursive algorithm for the coefficients in the 
"level" expansion (jA.Tjl 



i^o(A) = 1 (A.13) 

A- A, 



Hl{/\) = ^ . _"?" HL-mn{^m,n + mn) 



mn<L ' 



where we have again suppressed the dependence on Aj. 

Now for the purposes of quantum gravity we want to combine two blocks of different 
conformal field theories, conventionally be the "matter" and the "Liouville" one, with central 
charges respectively cu and cl- We do not necessarily require the "complementarity" of these 
quantities, introducing the "spectator" central charge 

Csp = 26-cm-cl (A. 14) 

to take care of the deficit. On the contrary we do require the complementarity of the "matter" 
and "Liouville" external dimensions Aj and Aj 

Ai + Ai = l (A.15) 

for i = 1,2,3,4. The "matter" and "Liouville" blocks are combined to 



(M) f Ai A3 
^ V ^2 A4 


\ _^(L) / Ai A3 
) A \^ A2 A4 


) 


xil-x)9l{q) 


(A.16) 


where we conventionally denoted 






^sp(g) = r"^/^(g^) 




(A.17) 


with 

CO CX) 








vi^) = 1 1 (1 


-<!'')= >;(- 


\n n{3n+l) 


(A. 18) 



k=l n=— 00 

the Dedekind function. 

Notice, that we do not demand the intermediate dimensions A and A to be complemen- 
tary. To avoid misunderstanding, let us stress that Hsp{q) is simply a convenient notation 
and hardly can be interpreted as a "contribution of the spectator matter" to the block and 
the four point function. Some additional simplifications of the product of the elliptic blocks, 
which occur if the "matter" and "Liouville" CFT's are indeed complementary (i.e., Cgp = 0) 
will be discussed in 
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